DENSITY FUNCTIONAL IN LIQUID HELIUM: RECENT 

> ■ APPLICATIONS TO DROPLETS AND SURFACES 

O 

. Franco Dalfovo, Andrea Lastri, and Sandro Strineari 

CN ' Dipartimento di Fisica, Universita di Trento 

1-38050 Povo, Trento, Italy 

0^ 
oo 
O 









1. INTRODUCTION 

In the present contribution we present some of the most recent results of density 
functional theory for inhomogeneous systems of liquid "^He. Almost ten years ago, 
one of us (S.S.) introduced a simple density functional for liquid *He and ^He. The 
functional was used to study the ground state properties of the free surface [1] and 



O ! of droplets [2]. The theory has been later refined and extended in order to improve 

^ ' its accuracy and to make it applicable to a more microscopic level. Finite range 

correlations have been explicitly introduced [3] to explore strongly inhomogeneous 
^^ . systems, like layered structures of liquid helium on solid substrates [4] and near 

H I impurities [5] or vortices [6]. ^He-^He mixtures have been also extensively studied 

with the same formalism [7]. 

It is a good time to see what the situation is. There are two major questions 
which have to be addressed. First, what is the present understanding of the foun- 
dations of density functional theory. For instance, what is the connection between 
phenomenological density functionals and microscopic approaches, as well as the ap- 
plicability of the theory to dynamic problems through a time dependent formulation. 
Second, what is the predictive power of the theory and its accuracy. In the present 
contribution we will mainly concentrate on the second problem. 

We will present a new density functional [8] (hereafter called Orsay- Trento func- 
tional), which is the result of a natural evolution and improvement of the theory 
starting from the original zero-range functional of Refs. [1,2]. We will focus on the 
key ingredients of the theory and we will show how they determine important pre- 
dictions concerning inhomogeneous systems. These predictions turn out to be quite 
accurate even at the level of interatomic distances. We will present a selected set of 
results, referring to Refs. [8,9] for a more detailed discussion. 



2. DENSITY FUNCTIONAL: EQUILIBRIUM PROPERTIES 

Let us consider the complex wave function 

vl/(r,t)=$(r,t)exp('^5(r,t)^ , (1) 

where the real function $ is related to the one particle density by p = $^, while 
the phase S is related to the velocity of the fluid by v = (l/m)V5'. In the density 
functional approach, the energy of the system at zero temperature is assumed to be 
a functional of \1/: 

E = fdrn[^,'^*] . (2) 



In the calculation of the ground state, only states with zero velocity are considered, 
so that the energy is simply a functional of the one-body density /o(r). The ground 
state is then found by minimizing the energy with respect to p. One usually write 
the energy as the sum of the kinetic energy of non-interacting particles ("quantum 
pressure" ) plus a term incorporating the effects of quantum fluctuations and corre- 
lations: 

E = jdrHoip] = /^r |--(Vv/p)2 + E,[p] . (3) 

Minimization with respect to the density leads to the Hartree equation 

{-|^V2 + t/[p,r]|v/;^ = PA^W , (4) 

where U[p, r] = dEc/6p{r) acts as a mean field, while the chemical potential p is 
introduced in order to ensure the normalization of the density. 

A derivation of functional Tio starting from first principle is not available. One 
then resorts to approximate schemes for the correlation energy. The phenomenolog- 
ical density functional approach consists of choosing a reasonable ansatz for Hq, in 
which a few parameters are fixed to reproduce known properties of the liquid. For 
example, the functional of Refs. [1,2] was written in the form 



E ^ dr 



£(VVP)^ + ^P^ + fp-- + <i(Vp)^ 



(5) 



where 6, c, 7 and d are phenomenological parameters fixed to reproduce the ground 
state energy, density and compressibility in the homogeneous liquid at zero pressure, 
as well as the free surface energy. Non locality effects have been included in a more 
adequate way by Dupont-Roc et al. [3], who generalized Eq. (5) accounting for the 
finite range of the atom-atom interaction and the existence of the hard-core. In that 
case, the term quadratic in the density is replaced by a two-body potential energy, 
containing a Lennard- Jones interaction, and the term with the power 7 is taken to 
depend on a "coarse grained density" , i.e., the density averaged on a sphere of atomic 
size. The three parameters of the functional are fixed to reproduce the bulk equation 
of state and the energy per particle at zero pressure. The most important feature 



Figure 1. Static response function in liquid "^He at zero pressure. Points: experi- 
mental data [10]; dotted line: from functional of Refs. [1,2]; dashed line: from Ref. [3]; 
solid line: from the Orsay-Trento functional. 



of this approach is that the static response function of the liquid turns out to be 
strongly g-dependent, with a peak at the roton wave length, in qualitative agreement 
with the experimental data [10]. 

The static response function x{q) plays a key role in density functional theory, 
since it fixes the response of the system to static density perturbations. It is conse- 
quently a quantity which, in principle, must be exactly accounted for by the theory. 
It is easily calculated from the density functional by taking the second functional 
derivative of the energy in g-space. Actually the functional of Ref. [3] provides a 
static response function which underestimates significantly the height of the peak in 
the roton region. For this reason we have corrected that functional, introducing a 
new non local term depending on gradients of the density. The new term is fixed to 
reproduce the experimental static response function in the liquid. The function Tio, 
entering Eq. (3), is then written in this way: 

Wo= ^(^Vpf + Ijdv' p{r)Vd\v-v'\}p{v') + |p(r)fe)2 + |p(r)fe)= 

-|la,/dr'F(|r-r'|)(l-^)vpW.VpM(l-^) . (6) 

The two-body interaction Vi is the Lennard- Jones interatomic potential, with the 
standard parameters a — 2.556 A and e = 10.22 K, screened at short distance 



(V = for r < hj with h = 2.1903A). The weighted density p is the average of 
p{r) over a sphere of radius h. Thus the two terms containing p, with the parameters 
C2 == -2.411857x10^ K A^ and C3 = 1.858496x10^ K A^, account phenomenologicaUy 
for short range correlations. The last term corresponds to a non local correction to 
the kinetic energy. The function F is a simple gaussian F{r) = 7r~^''^i~^ exp(— r^/£^) 
with e = 1 k, while a^ = 54.31 A^. The factors (1 - p/pos), with pos = 0.04 A-^, 
are included in order to obtain a pressure dependence of the static response function 
close to the one predicted by Diffusion Monte Carlo simulations [11]. In conclusion, 
the function TCq is such that the theory reproduces two important features of the 
liquid, namely the equation of state and the static response function. The situation, 
as concerns the quantity x{q) is shown in Fig. 1, where the experimental data are 
compared with the results of three different functionals. The q = limit is fixed 
by the compressibility of the system, which, in all cases, is an input of the theory. 
It ensures the correct behaviour in the long wavelength limit and, consequently, the 
correct description of systems characterized by slow density variations. The peak of 
the static response function in the roton region, g ^ 2 A~^, is very significant in 
characterizing structural properties on the interatomic length scale. 

With the functional described above we have studied the properties of the free 
surface, droplets and films. Here we restrict ourselves to a few relevant results for 
helium droplets, where a detailed comparison between the predictions of different 
density functionals with the ones of Monte Carlo calculations, is possible. In this 
case the Hartree equation (4) has to be solved in spherical geometry. The ground 
state energy and density of small droplets are shown in Fig. 2 and Fig. 3, respectively, 
as a function of the number of atoms. The density profile of a droplet of 70 atoms is 
shown in Fig. 4, together with the predictions of other calculations. Some relevant 
comments follow: 

• The predictions of density functional theory for the energy of droplets become 
closer and closer to the results of ab initio simulations, as the accuracy of the 
static response function is improved. We note that the three functionals give 
practically the same surface energy, very close to the experimental one, and the 
same bulk equation of state, while they differ in the behaviour of the static 
response function. The close agreement between the predictions of the Orsay- 
Trento functional and the ones of ab initio simulations seems to indicate that the 
new non-local kinetic energy term in Eq. (6) accounts properly for correlations 
on the atomic length scale. The detailed form of the functional is not crucial; 
different parametrizations give very similar results provided the static response 
function in bulk liquid is reproduced. 

• The density profiles exhibit small regular oscillations. These oscillations were not 
found with previous functionals, because the latter underestimated the height of 
x{q) at the roton wave length. The connection between the oscillations of the 
density profile and the behaviour of the static response function was suggested 
twenty years ago by Regge and Rasetti [14] . Subsequent microscopic calculations 
did not provide any conclusive evidence of this effect. The results of the density 
functional approach support the Regge's idea. These oscillations have a regular 
behaviour when one goes from small to very large clusters, up to limit of a 
planar free surface. The density oscillations are found to affect also the the 
evaporation energy, [E{N — 1) — E{N)], characterizing the mass distribution of 



Figure 2. Energy per particle of iielium droplets. Tiie results of the present Orsay- 
Trento functional (solid line) are compared with the ones of functional of Refs. [1,2] 
(dotted line) and of Ref. [3] (dashed line), as well as with Monte Carlo simulations 
of Ref. [12] (crosses) and Ref. [13] (circles). 



Figure 3. Density profile of helium droplets with A^ = 8 to 60. The density is 
normalized to the one of the bulk liquid. 



Figure 4. Density profile of a droplet witfi 70 "^He atoms. Solid line: present work; 
triangle: DMC simulations of Ref. [15]; dashed line: Variational (HNC) calculations 
[16]; circles: most recent DMC calculations by Barnett and Whaley [17]. 



droplets in experimental beams: deviations of the order of 0.1 K with respect 
to a pure liquid drop model are predicted. These deviations are signatures of a 
shell structure ("soft sphere close packing") induced by quantum correlation at 
the liquid surface. 
• Density oscillations have been recently found in Diffusion Monte Carlo (DMC) 
calculations on small droplets by Chin and Krotscheck [15], even if their results 
are probably affected by large statistical errors. The same authors have more 
recently employed a variational (HNC) method, and found oscillations closer to 
the ours [16]. A very good agreement with the density functional results has 
been very recently obtained by Barnett and Whaley in their most recent DMC 
calculations [17]. 

In other systems, like the free surface, films, porous media, doped clusters, less 
information are available from ab initio methods. The accuracy of the density func- 
tional approach is expected to be similar to the one shown for helium droplets. We 
emphasize again that the freedom in the choice of the parametrization of the start- 
ing functional is only a minor source of indetermination in the predictive power of 
the theory, provided the crucial ingredients (equation of state and static response 
function in bulk liquid) are properly accounted for. 



Table 1 

Values of the parameters used in Vj{r). 



dt I dr 



3. DENSITY FUNCTIONAL: EXCITED STATES 

Starting again from functional (2) one can develop a time dependent density 

functional theory in order to study the dynamics of inhomogeneous systems. In this 

case one has to apply the least action principle 

"^a r r 5^- 

n[^\-^]- ^i-^*-^ --^*ih— = . (7) 

The equations of motion for the excited states of the fluid can be derived by making 
variations with respect to \E' or \E'*. One finds a Schrodinger-like equation of the form 

(^-^)vE' = zn^vl> , (8) 

where H = dE/d"^* is an effective Hamiltonian. If one looks for linearized solutions 

^{r,t)=^o{r)+S^{r,t) (9) 

the Hamiltonian H takes the form H = Hq + 5H, where the static Hamiltonian Hq 
fixes the equilibrium state as in Eq. (4). The term 5H is linear in 5\1/ and accounts 
for changes in the Hamiltonian induced by the collective motion of the system. The 
Schrodinger equation (8) has to be solved using a self-consistent procedure. This 
linearized theory concides with the Random Phase Approximation for bosons. More 
detailed aspects of the theory were discussed by S. Stringari in the same Workshop. 
To solve the equations of motion of the time dependent density functional theory 
one has to consider explicitly the velocity dependence of the function Ti. We write 

n = nM+'>^v[pM , (lo) 

where Tio is the velocity independent function given in Eq. (6), while 7i„ is taken to 

be [8] 

^. = yP(r)|v(r)|2-^|rfr'K;(|r-r'|)p(r)p(rO [v(r) - v(r')]' . (11) 

The first term is the usual kinetic energy for free particles. The second term plays 
the role of a non local kinetic energy, accounting for backfiow effects through an 
effective current- current interaction Vj. A similar functional was introduced long 
time ago by Thouless [18] to study the fiow of a dense superfiuid. We fix Vj{r) using 
the phono n-roton dispersion in bulk liquid as input. We have chosen the simple 
parametrization 

yj{r) = (711 + 712^^) exp(-Q;ir^) + (721 + 722^^) exp(-a2r^) , (12) 

where the parameters are given in Table I. With this choice the phonon-roton spec- 
trum of the uniform liquid at zero pressure is well reproduced, and its pressure de- 
pendence also agrees with the experimental data. The theory is then used to predict 
dynamic properties of non uniform systems. 



Figure 5. Excitation spectrum of the free surface (see text) 



We show here what happens in the case of the free surface. The time dependent 
part of the wave function (9) has to be expanded in plane waves in the direction 
paraUel to the surface. The equations of motion are solved for the 2;-dependent 
coefficients, where z is the orthogonal coordinate. In practice, we expand the solutions 
on a basis of eigenfunctions of the static Hamiltonian, in order to write the equations 
in the form of a matrix diagonalization. One gets the transition density, the energy 
and the strength associated with each excited state, for any given value of the parallel 
wave vector k. The same procedure can be easily applied to the case of helium films, 
by adding the external potential of the substrate. 

In Fig. 5 we show the excitation spectrum of the free surface. Up to k ~ 1.15 
A~^, the state with lowest energy is localized at the surface. In the long wave length 
limit /c — !► its dispersion coincides with the hydrodynamic dispersion of ripplons 
(dotted line) . It deviates significantly from the hydrodynamic law when approaching 
the threshold energy A = 8.7 K for rotons. A similar deviation has been found 
in the dispersion of surface modes in helium films measured in neutron scattering 
experiments [19] (points with error bars). The bulk phonon-roton dispersion (upper 
solid line) is also plotted, together with the experimental data. We stress again that, 
while the phonon-roton dispersion in the bulk is used as input to fix the velocity 
dependent part of the functional, the theory becomes predictive for the dynamic 
properties of inhomogeneous systems, like the free surface. 

Below A the surface mode is undamped, while above A it couples with the con- 
tinuum of bulk modes (rotons with negative and positive group velocity) propagating 



at different angles {q^ 7^ 0). This results in a spreading of the strength associated 
with the surface mode. Actually the spreading predicted by our theory is rather 
small. This can be seen by plotting the dynamic structure function S{k,uj), defined 
by 



b 



^^'k. 



6{uj-uj]^^b) ■ (13) 



In Fig. 6 (top) we show the dynamic structure function for scattering at grazing angle 
{Qz = 0) on a slab (liquid between two parallel surfaces at a distance of 50 A). For 
graphical convenience we have replaced the 5-functions with gaussians of width 0.4 
K. The strength of the surface mode is well localized also above A, even though the 
strength is partially distributed among bulk modes coupled to ripplons. The position 
of the peak above A is shown also in Fig. 5 as a dashed line. The hybridization 
mechanism between ripplons and rotons is discussed on a general ground in Ref. [9] . 

The form of the spectrum can change significantly in liquid helium films on solid 
substrates, where the substrate-helium potential produces a layering of the helium 
density and the bulk modes are discretized by the finite size of the system along z. An 
example is given in Fig. 6 (bottom), where we show the dynamic structure function 
for a film on Sodium. The areal density of helium is 0.23 A~^, which corresponds 
to approximately 4 liquid layers. The surface mode is clearly visible. Near the 
roton minimum one notes a new structure coming from the first liquid layer close to 
the substrate. It is interpreted as a 2D roton, whose energy becomes significantly 
lower than A in the case of more attractive substrates, refiecting the tendency to 
solidification. Similar structures have been also seen in the experimental spectrum 
[19]. 

Finally, we note that the density functional method is also suitable to investi- 
gating the properties of refiection and evaporation of bulk excitations at the surface. 
Liquid helium is very peculiar from this viewpoint. For instance, one can study the 
one-to-one process of a roton impinging the surface and ejecting an atom (quantum 
evaporation). Experimental data about evaporation and condensation rates are now 
becoming available [20]. The threshold for quantum evaporation at zero tempera- 
ture is shown in Fig. 5 (dot-dashed line). The q = value is the chemical potential 
\n\ = 7.15 K. Above this curve the solutions of the time dependent density func- 
tional theory are combinations of free atom states, outside the liquid, and bulk states 
(phonons and rotons). From their coupling one can extract the evaporation and 
condensation probabilities. This work is in progress. 



4. CONCLUSIONS 

In this contribution we have briefiy traced the evolution of phenomenological 
density functional theories, from the one of Refs. [1,2] to the most recent Orsay- 
Trento functional [8]. We have emphasized the crucial role of the static response 
function of the bulk liquid; together with the equation of state, this quantity has 
to be considered the main ingredient of the density functional apprach. We have 
shown how the predictions of the theory improve when the static response function is 
properly reproduced. For the time dependent version of the theory one can introduce 
backflow effects through an effective current-current interaction. In that case the 



Figure 6. Dynamic structure function for a 50 A thick slab (top) and a film on 
Sodium (bottom), with coverage of 0.23 A~^. 



phonon-roton dispersion in the uniform liquid is used as input. This makes the 
theory quantitative in the prediction of the excited states of inhomogeneous systems. 
We have shown results for the dispersion of ripplons on the free surface and for the 
excitations of films. An accurate description of the excitations in bulk is particularly 
important to study the ripplon-roton hybridization mechanism [9] and the process of 
quantum evaporation. 

The predictions of the theory compare well with available results of Monte Carlo 
similations in small systems. The density functional approach can be easily applied 
to large systems too, and consequently can be considered a reliable tool to exploring 
the properties of inhomogeneous liquid helium in a rather systematic way. 
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